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ARTICLE INFO ABSTRACT 
Article history: The gastrointestinal peptide hormone ghrelin promotes food intake and increases body weight and 
Received 15 June 2010 adiposity through activation of the growth hormone secretagogue receptor (GHSR1a). To promote its 
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biological action ghrelin is acylated at its serine 3 residue by the recently discovered ghrelin O-acyl- 
transferase (GOAT, a.k.a. membrane-bound O-acyltransferase 4, MBOAT4). Plasma levels of total and acyl- 
ghrelin are negatively correlated with body-mass-index (BMI); as lower the BMI as higher plasma levels 
of total and acylated ghrelin and vice versa. Accordingly, plasma levels of total and acyl-ghrelin are 


Silas elevated in patients with anorexia nervosa (AN) and decline upon weight regain. The importance of the 
Ghrelin endogenous Goat/ghrelin system in the neuroendocrine adaptation to fasting was recently highlighted 
Anorexia nervosa by the observation that acyl-ghrelin mediated elevation of growth hormone (GH) release prevents 
Association study starvation induced hypoglycemia in Goat~!~ mice. The aim of this study was to test if genetic variation of 
GOAT is implicated in the etiology of AN. We therefore assessed association of 6 tagging single nucleotide 
polymorphisms (tagSNPs), which were predicted to cover 96% the common genetic variability of GOAT 
plus 50 kb of the 5’ and 3’ flanking region, in 543 German patients with AN and 612 German normal and 
underweight healthy controls. Based on a recessive mode of inheritance we observed some evidence for 
association of the G/G genotype at SNP rs10096097 with AN (nominal two-sided p = 0.031). Based on our 

results we conclude that genetic variation in GOAT might be implicated in the etiology of AN. 
© 2010 Elsevier Ltd. All rights reserved. 

1. Introduction 

Anorexia nervosa (AN) is an eating disorder of unknown 
~* Corresponding author. Tel.: +-1 513 558 0433. etiology, affecting approximately 1—2% of adolescents, mostly 
E-mail address: tschoemh@ucmail.uc.edu (M.H. Tschop). young women between the ages of 14 and 18 years (Vyver et al., 
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2008; Hoek and van Hoeken, 2003). Highly efficacious or even 


curative pharmacological treatment options remain to be discov- GSSrALTSS Sr ery sRgesena 
A : adi . AAGoHToRMnDAWeR Sr OAN 
ered. Even though substantial evidence indicates that approxi- SSSseoSdSsSsSSdcSSdSSG 
mately 50—80% of the variance of AN can be explained by genetic Se ieee eee EN OO ee 
factors (Klump et al., 2001), little is known about the molecular 
mechanisms of this disorder. Several studies have suggested that y 
genes implicated in the regulation of food intake and energy sg mS 
metabolism might be implicated in the etiology of AN (Hebebrand & oi ra ea iA be 5 
and Remschmidt, 1995; Hinney et al., 2004; Pinheiro et al., 2006; = ba 7 i 8 a 8 = 
Sulek et al., 2007; Scherag et al., 2010). . 2 
The gastrointestinal peptide hormone ghrelin promotes body 3 E 
weight gain and adiposity by leading to increased food intake and § ee) 88 22 BS BS RR BIE 
decreased fat utilization in both rodents (Tschép et al., 2000; > oI x oA m 3 6 m % = 3 S § Q + 
Asakawa et al., 2001) and humans (Wren et al., 2001; Hotta et al., 3 meee a Os ee a ae 3 
2009; Druce et al., 2006). Numerous studies are currently revealing a 5 
ghrelin’s mechanisms of action, which are believed to rely on specific 3) Ree a on ooORe ene oe one = 
hypothalamic as well as extra-hypothalamic neurocircuitries. Until 2 = a i = 3 iB ae i a 3 wu 3 ae S = 
recently however, relatively little was known about the mechanisms 3 a . iB - = ef = _ = = - S. = = = = = ~ S 
of ghrelin activation, as the hormone’s ability to bind and activate its = BAPE? RAPS Fo RAY rs 
known receptor (GHSR1a) depends on post-translational acylation of a i) 
the ghrelin peptide at the serine 3 residue (Kojima et al., 1999). e = g 
Recently, the membrane-bound O-acyltransferase 4 (MBOAT4), 8 = a 
which was later renamed to ghrelin O-acyltransferase (GOAT), was gS g 5 
identified as the endogenous key enzyme responsible for ghrelin E 2 4 
octanoylation (Yang et al., 2008; Gutierrez et al., 2008). The obser- io THOUUCUUKEUUECUUECUUUYUY|S 
: . eave : N S229 SZC SZ SG E2GS29NUE|F 
vation that acylated ghrelin is completely absent in mice lacking Goat S|, ae 
(Goat!) indicates that Goat is the only enzyme capable of acylating Ee 2 a ; 
ghrelin in-vivo (Kirchner et al., 2009). Sisl|_e = = ae = = ox 
f . r ze}OlZo © © © © © Bs 
In rodents, plasma levels of ghrelin are increased in response to = a 
fasting and decrease upon re-feeding (Tsch6p et al., 2000; Asakawa 8 ae 
et al., 2001; Kirchner et al., 2009). Accordingly, plasma levels of gy 3 2 5 
total and acyl-ghrelin are elevated in patients with AN (Otto et al., g 4 ia 5 
2001; Nakahara et al., 2008; Germain et al., 2007, 2009; Prince 5 a9 S © x N > +8 
et al., 2009; Monteleone et al., 2008; Tolle et al., 2003) and a = Ss 3 S a S = £3 
decline to normal values upon weight restoration (Otto et al., 2001, ¢ 8 5 
2005; Troisi et al., 2005). Interestingly, several small studies 2 sees Gi, Gig bse oA ages, elle e 
reported higher plasma levels of acyl-ghrelin in patients with AN g Oo ace Re of 2a eo 3 5 g 
when compared to BMI-matched lean controls (Germain et al., a a ee eee ee ee, Gee oy es 
2007, 2009; Tolle et al., 2003). S| ee ee ee 
Ghrelin promotes food intake through stimulation of hypotha- o % 5 
lamic orexigenic neuropeptides, such as neuropeptide Y (Npy) and rf oe ae ae E 8 
agouti-related protein (Agrp) (Kamegai et al., 2000, 2001). Despite i SRASSSSFASERSSISESSISE 
its orexigenic effects ghrelin promotes the expression and release of ‘y ZYVSQLVSHVengatrreterealy 3 
growth hormone (GH) from the anterior pituitary gland and thus ¢, on eT a ae “ . z 
indirectly triggers the release of hepatic insulin-like growth factor-1 < 2 A 2 
(IGF-1) (Kojima et al., 1999; Arvat et al., 2001 ). Both, GH and IGF-1 are e 3g SYYSBVSBVsUIsUVyyy = 3 
anabolic hormones that increase lean body mass by stimulation of = 3 3 Bo 
skeletal muscle growth and inhibition of skeletal muscle protein S| 8 S35 
breakdown (Fuller et al., 1992; Moller et al., 2009). Resistance to GH Zz 2/22 + * 2 + + s ro 
is typically accompanied by high plasma levels of GH and low plasma 2 ae 
levels of IGF-1 (Hanaire-Broutin et al., 1996; Monzavi and Cohen, z | & 
2002). Impairment of the GH/IGF-1 axis is frequently described in = § ee ee a) ee ee ee Py g 
patients with AN; plasma levels of GH are typically elevated while a a) > S > 2 > E a 
IGF-1 levels are decreased (Counts et al., 1992; Scacchi et al., 1997; a ae _ = a ae ” a2 
Misra et al., 2003). The importance of the Goat/ghrelin system in = ee 
the modulation of the GH/IGF-1 axis was recently demonstrated by 2 9g 9g C S a i 2 € 
Zhao and colleagues. They showed in Goat~!~ mice that ghrelin 2 7 = fe <x < < oO z Pi 
mediated elevation of GH prevents death by preserving blood 3. ge a = 9 8 = o8 
glucose levels under conditions of severe caloric restriction (Zhao & =z 7 N Ss Ss ba g g 
et al., 2010). Together, the observation that plasma levels of total a DS 
and acyl-ghrelin are increased and that the GH/IGF-1 axis is s * om 5 3 
impaired in patients with AN indicates that genetic variation of GOAT a 3 8 = > B 2 3 ‘y 
might be implicated in the etiology of AN. The aim of the present = 8 g a] & ip a ra Es 
case-control association study was therefore to assess the potential Ze Z 3 4 7 2 : 2 Za 
implication of genetic variation in GOAT with AN. ES 8 
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2. Methods 
2.1. Study subjects 


The ascertainment strategy was previously described in detail 
(Hinney et al., 1997, 2006). Written informed consent was given by 
all participants and in the case of minors, by their parents. The 
study was approved by the Ethics Committees of the Universities of 
Marburg and Duisburg-Essen, Aachen, Tiibingen, Heidelberg, 
Erlangen, Freiburg, Wiirzburg, Bochum and the Charité Berlin, and 
carried out according to the Declaration of Helsinki. 

Study group 1 (AN cases) comprised 543 (1 male) German 
patients with AN. All patients with AN were recruited based on the 
DSM-IV criteria for anorexia nervosa (American Psychiatric Asso- 
ciation, APA, 1994). Of these, 460 (1 male) individuals had acute AN, 
whereas 81 (0 males) were recovered individuals from long-term 
follow-up studies of AN (Herpertz-Dahlmann et al., 1996, 2001; 
Zipfel et al., 2000). At inclusion, the acute patients had a mean 
age of 19.79 + 7.80 years and a mean BMI of 15.68 + 1.74 kg/m. The 
recovered patients had a mean age of 27.89 + 9.49 years and a mean 
BMI of 20.72 + 2.57 kg/m”. In total, the 543 patients with AN had 
a mean age of 21.00 + 8.57 years and a mean BMI of 16.44 + 2.61 kg/ 
m7 (mean BMI percentile based on the German National Nutrition 
Survey (Hebebrand et al., 1994; Anders et al., 1990) of acute and 
recovered AN patients was 3.49 + 7.93 and 38.80 + 28.78, 
respectively). 

Study group 2 (controls) comprised 612 German normal and 
underweight healthy controls (229 males) with a mean age of 
25.66 + 5.13 years and a mean BMI of 19.11 +1.91 kg/m? (mean BMI 
percentile 15.87 + 19.55) (Miiller et al., 2010). All controls with 
a current eating disorder like AN or Eating Disorder not Otherwise 
Specified (EDNOS), as assessed by a short version of the composite 
international diagnostic interview (CIDI) were excluded from the 
study. The CIDI questionnaire was kindly provided by Prof. 
Wittchen. 


2.2. Molecular genetic methods 


Using the ‘Algorithm-Tagger-Pairwise Tagging’ approach 
implemented in HapMap (www.hapmap.org, HapMap Data Release 
24/Phasell Nov08, NCBI36 assembly, r? cut-off: 0.8, minor allele 
frequency (MAF): 0.1), six tagSNPs (rs1355412, rs10096097, 
1816876504, rs3735989, rs13272159 and rs4733400) were pre- 
dicted to cover 96% of the common genetic variability of the GOAT 
coding regio + 50 kb of the 5’ and 3’ flanking region as derived from 
the HapMap-CEU population. The analyzed tagSNPs are located in 


Table 2 
Primer sequences and restriction enzymes. 
SNP Primer (5’ - 3’) 
rs10096097 F: 5’-AGC CAG TAA GGC AGT TCC AA-3’ 


R: 5’-TCA CCA CGG AAA ACA TTG AA-3’ 


1s4733400 F: 5/-GGT GGA GAG TCA GCT TTG CT-3’ 
R: 5’-ACA CAG CAT GAT CCC GTC TC-3’ 
rs1355412 F: 5’-GAG CAA TGG GCA GAG TCT TG-3’ 
R: 5’-GAT GGC CAG GGG TAA TIT CT-3’ 
rs13272159" F-in: 5’-AGG AAA CAT TIT AGA AAG TGC CAG ATT AC-3’ 
R-in: 5’-AAA GTT CAA ATG TAT GCA TTA TCG ATT CA-3’ 
F-out:5’-AGT GTC TIT CCC ATG GTG TAC TGT ATT C-3’ 
R-out:5’-GTA AAC CTT AAG GGA AGA AAC CAC ATC T-3’ 
rs16876504 F: 5’-AGA GCA TCA GCC AGG ACATT-3’ 
R: 5’-AAC CAC CTA ATA CAG TCA CAG TTG A-3’ 
183735989 F: 5’-ACG CTT GCT AGC AGC ATT TT-3’ 


R: 5’-GGA GGC AGA GGT TGC AGT AA-3’ 


the 5’ flanking region (rs10096097, 1rs4733400), exon 3 
(rs13272159) or the 3’ flanking region (rs1355412, rs16876504, 
1s3735989) of GOAT (Table 1). 

Genotyping was performed using restriction-fragment length 
polymorphism analyses (PCR-RFLP). Only SNP rs13272159 was 
genotyped using tetra-ARMS-PCR (Ye et al., 2001) as described 
previously (Miiller et al. 2008a). The corresponding primer 
sequences and restriction enzymes are displayed in Table 2. For 
validity of the genotypes, alleles were rated independently by at 
least two experienced individuals. Discrepancies were resolved 
unambiguously either by reaching consensus or by retyping. 
Additionally, several random samples of patients being homozy- 
gous for the low frequency alleles of SNPs rs10096097 and 
1813272159 were resequenced to confirm reliability of genotyping. 


2.3. Statistics 


The genotype distribution of the analyzed variants was tested 
for deviations from Hardy—Weinberg equilibrium using an exact 
two-sided test (Wigginton et al., 2005). As we had no specific 
a priori hypothesis for the underlying mode of inheritance we used 
the Cochran—Armitage trend test for the log-additive, recessive and 
the dominant mode of inheritance (Armitage, 1955) to test for 
association with AN. 

Additionally, odds ratios (OR) and 95% confidence intervals (95% 
Cl) were estimated based on a logistic regression with and without 
adjustment for sex for a 2-parameter model and three 1-parameter 
models for which each genotype was compared to the other 
genotypes. All reported p-values are two-sided, nominal, not 
adjusted for multiple testing and if not stated otherwise based on 
asymptotical considerations. 

Power calculations were performed using the program QUANTO 
version 1.2.4 under the assumption of a MAF of 0.3 as observed in 
the HapMap-CEU sample for rs10096097. 543 cases and 612 
controls were estimated to yield a power >0.80 to detect an OR of 
1.28, 1.68 and 1.39 for a log-additive, recessive and dominant model 
with the minor allele as risk allele, respectively. 


3. Results 


Analysis of 6 tagSNPs in 1155 individuals (543 German patients 
with AN and 612 German normal and underweight healthy 
controls) revealed some evidence for an association of the G/G 
genotype at SNP rs10096097 with AN (recessive model; nominal 
two-sided p = 0.031, Table 1). The corresponding genotypic odds 
ratio (OR) and 95% confidence interval (95% CI) for the G/G 


Product size Enzyme? RFLP-PCR Fragment size 
153bp Hhal G: 114bp + 39bp 
A: 153bp 
235bp Dral A: 133bp + 102bp 
G: 235bp 
396bp Hpy188I G: 145bp, 128bp, 96bp, 27bp 
A: 241bp, 128bp, 27bp 
150bp Faul G: 104bp + 46bp 
A:150bp 
232bp Btscl C: 187bp + 45bp 
G: 232bp 


* Product size of outer primer is 302bp, product size for C and T-allele is 203bp and 157bp, respectively. 


> all enzymes have been purchased from NEB, Germany. 
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Odds ratios and 95% confidence intervals for SNP rs10096097 of GOAT based on a logistic regression analysis with and without adjustment for sex. 


Table 3 
Genotype/Model OR 95% CI 
A/A 1 
A/G 0.788 0.616—1.009 
G/G 1.490 0.940—2.361 
Recessive Model: G/G vs A/G + A/A 1.636 1.043—2.566 


* Sex adjusted odds ratio (OR). 
> Sex adjusted 95% confidence interval (CI). 
© Sex adjusted p-value. 


genotype compared to G/A and A/A based on a logistic regression is 
OR 1.613, 95% CI 1.037—2.508, p = 0.032, Table 3. This effect was 
stable after stratification for sex (Table 3). Noteworthy, the associ- 
ation of the G/G genotype at SNP rs10096097 with AN was stable 
after exclusion of male subjects (data not shown). 

Notably, a trend towards deviation from Hardy—Weinberg 
equilibrium was observed for genotype frequencies of SNPs 
rs10096097 and rs13272159 in AN cases (p = 0.002 and p = 0.016, 
respectively) but not in controls (p = 0.43 and p = 0.48, Table 1). 
This finding, however, is not surprising and expected in case of 
a true genetic association and indeed the number of homozygotes 
for the rs10096097 risk G-allele was increased in AN patients 
compared to controls (frequency of risk G/G genotype in AN cases 
and controls is 9.24% and 5.94%, respectively; Table 1). However, to 
ensure validity of genotyping all homozygous G-allele carrier at 
rs10096097 and rs13272159 were re-genotyped. Additionally, 
several random samples were resequenced. The genotypes were in 
100% accordance. 


4. Discussion 


Ghrelin is the only circulating molecule to significantly promote 
food intake and body weight, which it achieves through activation 
of the ghrelin receptor (GHSR1a). Recently, the ghrelin O-acyl- 
transferase (GOAT) was identified as the endogenous key enzyme 
responsible for ghrelin activation (Yang et al., 2008; Gutierrez et al., 
2008). 

The human GOAT gene is located on chromosome 8p12 and 
comprises three exons spanning 13.02kb. The human transcript 
contains 1.308kb which are translated into a protein of 435 amino 
acids from which Asp307 and His338 are essential for the enzymes 
ability to octanoylate ghrelin (Yang et al., 2008; Gutierrez et al., 
2008). The GOAT protein is highly conserved between mammals 
and the expression is largely restricted to ghrelin expressing 
tissues, such as the stomach and the intestine (Yang et al., 2008; 
Kirchner et al., 2009). Due to its key role in mediating ghrelin 
acylation, GOAT may offer potential in the treatment of eating/ 
wasting disorders and cachexia. 

Several studies have shown an inverse relation between the 
body-mass-index (BMI) and circulating levels of acylated ghrelin 
(Tolle et al., 2003; Nakai et al., 2004; Harada et al., 2008). Accord- 
ingly, compared to normal weight healthy controls, plasma levels of 
acyl-ghrelin are elevated in patients with AN, even during all 
phases of an oral glucose tolerance test (Nakai et al., 2004; Harada 
et al., 2008). 

Here we describe to our knowledge, the first analysis of the 
potential relevance of genetic variation in the human GOAT gene 
and anorexia nervosa. Analysis of 6 tagging SNPs in 543 German 
patients with AN and 612 German normal and underweight healthy 
controls revealed some evidence for association of the G/G geno- 
type at SNP rs10096097 with AN (recessive model, nominal 
p = 0.031). The corresponding logistic regression analysis showed 
that homozygous risk G-allele carriers have an increased risk for 
AN. However, the potential implication of SNP rs10096097, which is 


p-value Adj. OR? Adj. 95% CIP Adj. p-value® 
1 

0.059 0.824 0.623—1.090 0.174 

0.090 2.014 1.121—3.319 0.019 

0.032 2.172 1.221—3.862 0.008 


located 24.981bp upstream the translation start codon (ATG) of 
GOAT, for AN remains unknown. Gain-of-function mutations in 
GOAT might, at least in part, account for the observed higher plasma 
levels of acyl-ghrelin and growth hormone in patients with AN. 
Persistent hyperghrelinemia might entail ghrelin and growth 
hormone resistance, comparable to insulin resistance in persis- 
tently hyperinsulinemic individuals. The observation that patients 
with AN show higher plasma levels of acyl-ghrelin compared to 
BMI-matched healthy controls (Germain et al., 2007, 2009; Troisi 
et al., 2005) supports this hypothesis. On the other hand loss-of- 
function mutations in GOAT might lead to impaired eating behavior 
and thus an increased risk for AN. 

Several association studies for AN have previously focused on 
genes implicated in the regulation of energy balance, such as on the 
genes coding for ghrelin (Monteleone et al., 2006; Cellini et al., 
2006; Dardennes et al., 2007; Gueorguiev et al., 2009), the ghre- 
lin receptor (Gueorguiev et al., 2009; Miyasaka et al., 2006), leptin 
(Hinney et al., 1998), the leptin receptor (Quinton et al., 2004), the 
brain derived neurotrophic factor (BDNF) (Ribasés et al., 2003, 
2004, 2005; Friedel et al., 2005), the tumor necrosis factor-alpha 
(Ando et al., 2001) or genes of the endocannabinoid system (Miiller 
et al., 2008b). Most of these studies have yielded negative results. 

Eating disorders are typically characterized by diagnostic 
instability, and migration between different subtypes, and even 
between different eating disorders (Milos et al., 2005; Tozzi et al., 
2005). Consistency of diagnostic criteria over several months is 
therefore a requirement to clearly distinguish between different AN 
subtypes. Accordingly, one of the limitations of our study is that 
detailed amnesic information, especially regarding previous eating 
disorders as well as follow-up data, were not available for all 
samples. However, whereas several studies consistently show that 
25-30% of patients with eating disorders migrate from AN to BN, 
migration from BN to AN is less frequently described (Tozzi et al., 
2005; Gorwood et al., 2003). Accordingly, it is unlikely that the 
observed association of the rs10096097 G/G genotype with AN can 
be explained by AN patients with a history of BN. In summary we 
here report some evidence for an association of the G/G genotype at 
SNP rs10096097 5’ to the GOAT gene with AN. Additional studies in 
larger samples are warranted to confirm our data. 
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